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We present spin-resolved measurements of the local den- 
sity of states in Si doped GaAs. Both spin components ex- 
hibit strong mesoscopic fluctuations. In the magnetic quan- 
tum limit, the main features of the spin-up and spin-down 
components of the local density of states are found to be iden- 
tical apart from Zeeman splitting. Based on this observation, 
we introduce a mesoscopic method to measure the jr-factor in 
a material where macroscopic methods are severely restricted 
by disorder. Differences between the spin-up and spin-down 
components are discussed in terms of spin relaxation due to 
spin-orbit coupling. 



Spin-polarized electronic transport is currently attract- 
ing a lot of interest from both a fundamental and an 
applied point of view 0| . Spin-polarized transport is es- 
sential for the operation of spin transistors or spin valves, 
and recently the spin injection from a ferromagnet into 
a semiconductor could be observed |^ . Spin effects have 
also been proposed to be utilized for quantum computa- 
tion. In this respect it is important to study spin relax- 
ation and coherence in semiconductors |3j. 

Particularly interesting spin phenomena were observed 
in mesoscopic semiconductor structures. Prominent ex- 
amples are the Kondo effect in quantum dots and 
spin-polarized tunneling through impurity levels. Reso- 
nant tunneling through impurities is also an established 
technique to image the local density of states (LDOS) of 
doped GaAs 0-0. The LDOS exhibits mesoscopic fluc- 
tuations which can be understood in terms of interference 
of elastically scattered electron waves However, a 

spin-resolved measurement was impossible to date. 

In this paper, we investigate the spin dependence of 
the LDOS via resonant tunneling through a spin-split 
impurity level. The strongest fluctuations of the spin- 
up and spin-down components of the LDOS are found to 
be identical apart from Zeeman splitting. This observa- 
tion allows us to mesoscopically determine the gi-factor 
in a material where macroscopic methods - such as pho- 
toluminescence or magnetotransport - are restricted by 
disorder. Differences between the spin-up and spin-down 
components are discussed in terms of spin relaxation due 
to spin-orbit coupling. 

Our experiment is based on a strongly asymmetric 
double-barrier heterostructure, which was grown on n~^- 
type GaAs substrate. It consists of a 10 nm wide GaAs 
quantum well and two Alo.3Gao.7As barriers of 5 and 8 
nm width. The nominally undoped active region is sand- 
wiched between 300 nm thick GaAs contact layers doped 



with Si. The donor concentration has been experimen- 
tally determined to 3.3 x 10^^ cm^^ From this ma- 
terial we fabricated a 2 fim diameter mesa with Ohmic 
contacts, which contains a small number of impurities in 
the quantum well. 

We use the energetically-lowest impurity state S as 
spectrometer for the LDOS i' in the emitter contact 
adjacent to the thick barrier, see Fig. |[ In a mag- 
netic field B, the spectrometer exhibits a spin splitting 
AEs = gsf^sB, with gs the 5-factor of the impurity 
and fiB the Bohr magneton. The spin splitting in the 
emitter is AE =^ gfisB. By applying a bias voltage 
V, the spin-split levels Si and Si are shifted in energy 
with respect to the emitter. They cross the Fermi level 
at Vs^ i = {Es-j I — iJi'^^)/ea and two current steps re- 
sult from resonant tunneling. Here, Es^ and Es^ are 
the energies of S-\ and Si, /i^ is the chemical poten- 
tial in the emitter (measured from the band edge in the 
well and emitter, respectively), and a is the voltage-to- 
energy conversion coefficient. These two current steps 
can be clearly observed in Fig. E](a), which shows the 
current-voltage characteristic I{V) of our device at the 
base temperature of T = 20 mK of a dilution refriger- 
ator. The spin splitting between the steps increases as 
B II / is increased from 8.0 to 14.2 T. The overaU shift of 
both steps to higher bias voltage is a consequence of the 
diamagnetic shift of Es to higher energy. 

As the tunneling rate of the thick emitter barrier is 
much lower than that of the collector barrier, the tun- 
neling current images the LDOS in the emitter below 
the Fermi level /i^ at the position of |0|. If spin 
is conserved during tunneling, the current on the short 
plateau of the first step is proportional to the spin-up 
component of the LDOS, Ji cx while the current 
on the plateau of the second step reflects a superposi- 
tion of both spin components, I2 vi + At a bias 
voltage of V|_| = {Es^ ^ — E)/ea, the spin- up and spin- 
down components of the LDOS at energy E are probed 
by our experiment (measured from the band edge in the 
emitter). The oscillatory fine structure on the current 
plateaus in Fig. ^a) reflects LDOS fluctuations in the 
heavily doped emitter |Q. These fluctuations are more 
pronounced in the numerically obtained differential con- 
ductance G — dl/dV shown in Fig. 0(b), where sharp 
peaks result from the current steps in the I{V) data. 

To determine the g-f actor of the spectrometer, we plot 
AVs — Vs^ — Vs^ as a function of magnetic held. Fig- 
ure H shows that the data follows a straight line, which 
confirms that the two steps in the I{V) curves refiect 
one spin-split impurity level. From the temperature de- 
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pendence of the current step at B = T (not shown), 
we deduced a vohage-to-energy conversion coefficient of 
a = 0.50 for our device mTl|,|lJ]. Using this value, we 
obtain \gs\ = ea{dAVs /dB)/ fiB = 0.14. The ^-factor 
of the impurity is different from the bulk value due to a 
change of the electron energy as a result of the double- 
barrier confinement and the penetration of the electron 
wavefunction into the barrier material [Q . Note that the 
spin splitting of the spectrometer does not tell whether 
gs is positive or negative. 

Now we turn our attention to the LDOS in the mag- 
netic quantum limit, in which only one spin-split Landau 
band remains populated in the emitter. The transition 
to this limit occurs at B ~ 11.4 T. Figure ||(a) shows 
a color-scale image of the current as a function of bias 
voltage and magnetic field. It is dominated by two sharp 
current changes denoted by Vs^ and V^^ , which reflect 
the spin-split steps in the I(V) curves in Fig. |l](a). The 
magnetic-field dependence of Vs^ has been eliminated by 
shifting the I{V) curves in voltage. Between Vs^ and Vs^ 
as well as above Vs^ , we observe pronounced fan-type 
resonances which reflect the local structure of Landau 
bands in the heavily doped emitter |^ . These resonances 
in the LDOS are formed by quasi-one-dimensional inter- 
ference of multiply-scattered electron waves subject to 
Landau quantization. The slope of the LDOS resonances 
is dYs/dB = 3.0 mV/T |l|l. 

To study the spin dependence of the LDOS, we define 
two cuts at Vi and V2 [dashed white lines in Fig. ||(a)]. 
Figure ^(b) shows the corresponding currents Ii and I2 
as a function of magnetic field. Assuming in a first at- 
tempt that the spin- up and spin-down components of the 
LDOS in the emitter are identical apart from a Zeeman 
shift AE, we can recalculate hiB) from Ii{B) as follows. 
From the sketch in Fig. |4| we deduce 



from Eq. (H). Using Eq. (|l|), we now rewrite Eq. (^) into 



and 



h{B)<xu^{Ei,B) (1) 



l2{B)(x:^^iE2,B) + :^l{E2 + AEs,B) (2) 



where Ei and E2 correspond to Vi and V2 and AEs 
denotes the spin splitting of the spectrometer. Using 
vi{E,B) — — AE,B), we replace the spin-down 

component of the LDOS in Eq. by its spin-up coun- 
terpart. In addition, we exploit the linear, fan-type char- 
acter of the LDOS resonances and use z/|(i5i + E,B) = 
v^{Ei^B + [dB /dEE]E) to rescale all energy values to 
El. Defining £'21 — E2 — Ei, we obtain 

dB 



l2{B)^u^ [Ei,B + -^E2i 



Ei,B + —-[E21 + AEs - AE] (3) 



l2{B)=jh(B 



dB 
dE^ 
dB 



E21 



ih [E21 + AEs -AE]j . (4) 

Here, 7 was introduced to compensate a difference in the 
coefficients of proportionality in Eqs.(^ and (^. It takes 
into account that the impurity can only be occupied by 
one electron at a time due to Coulomb interaction 
For our analysis, we estimate 7 as {I2) b /"^{h) b — 0.98 
where (...)_b symbolizes averaging over magnetic field. 
The slope dB/dEE = l/[ea{dVE I dB)] = 0.67 T/meV 
and \gs\ — 0.14 have been measured and, in a first step, 
we assume g = —0.44 ||l^. 

Figure ^ compares the experimental data for I2 with 
I2 recalculated according to Eq. (|^). Assuming in addi- 
tion gs = —0.14, the strongest features in the LDOS are 
nicely reproduced by our reconstruction. Thus we con- 
clude that the spin dependence of large fluctuations in the 
LDOS can be understood in terms of Zeeman splitting. 
These features vary slowly in magnetic field (on a scale of 
about 0.5 T or larger). The agreement also corroborates 
that both the spin-up and spin-down components of the 
LDOS are measured at the same location. The recon- 
struction does not work for 175 — -1-0.14, which strongly 
suggests that < as predicted by calculations of the 
,g-factor in a 10 nm wide quantum well ||lj,|l^ . Strong de- 
viations between I2 and I2 are observed for features that 
vary rapidly in magnetic field (on a scale shorter than 0.1 
T). This is more evident at the bottom of Fig. ^, where 
we plot A/2 — I2 ~ I2 vs magnetic field. 

In order to determine the actual g-factor in the heavily 
doped GaAs emitter, we minimize the root-mean-square 

deviation cta/ = {(AhiB) - {Al2)Bf)T with respect 
to g. Figure ^ shows a clear minimum at g ~ —0.35. 
The accuracy of g is limited by the uncertainty of gs- 
Assuming an error of ±0.01 in gs, we calculate the error 
of g to be smaller than ±0.02. In our experiment, the g- 
factor is measured at the Fermi energy of about 25 meV 
above the band edge and at magnetic fields of 10.5 — 
15 T, while the standard value g = —0.44 was obtained 
at the band edge and at zero field The observed 

reduction of \g\ is in accordance with theory predicting a 
reduction of the 5-factor in high magnetic fields and at 
finite electron densities |T7| ]. 

Now we address the origin of the difference between I2 
and I2 in Fig. ^. One issue to be considered is energy 
relaxation. The "holes" generated by tunneling in the 
emitter Fermi sea have an energy deficit ea{V — Vs^ ^ ) 
with respect to he (see sketch in Fig. |l|). The higher the 
bias voltage the larger is this energy deficit, resulting in 
a faster relaxation via electron-electron scattering. In- 
elastic scattering reduces the phase-coherence length in 
the emitter and thus suppresses the LDOS fluctuations. 
This effect is clearly visible in Fig. |lj where the ffuctua- 
tions become less pronounced with increasing bias volt- 
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age. Similarly one might expect that I2 is different from 
I^, since I2 is measured deeper below the Fermi level 
than Ii which is used to calculate /|*. To check this idea, 
we calculated A/2 — I2—I2 f^'" 0-5 mV < V2 < 1 mV cor- 
responding to an energy range of 0.25 meV. Interestingly, 
A/2 is fairly independent from V2, which demonstrates 
that energy relaxation is irrelevant in our experiment. 

Thus the difference between I2 and I2 either indicates 
an inaccuracy in our reconstruction - possibly resulting 
from a magnetic-field dependence of the parameters 7, 
dB/dEE, gs, g - or it gives evidence of additional spin 
effects beyond Zeeman splitting. While we will not be 
able to give a definite answer, it is worth while to examine 
the influence of extra spin effects on the LDOS. 

One possibility is spin relaxation due to spin-orbit cou- 
pling in the emitter H. In this case, we can estimate the 
spin-relaxation time as follows. Semiclassically speaking, 
the LDOS fluctuations are a result of quantum interfer- 
ence of diffusive electron waves at length scales between 
the mean free path and the phase-coherence length 0. 
Long interference paths correspond to weak features that 
vary rapidly in magnetic field. By calculating the dif- 
ferential conductance G = dl/dV, we favor small and 
rapid variations resulting from long interference paths 
[see Fig. |l|(b)]. Due to the limited energy resolution hT 
of our spectrometer, the longest and dominant length 
scale in the conductance is = ■^/-D/r, where D is the 
diffusion constant in the emitter The correlation 

magnetic field of the conductance fluctuations is 

B,,G-^o/ll^4'o/{D/T) . (5) 

By calculating the difference A/2 = I2 — I2 , we elimi- 
nate the spin-independent features of the LDOS. Hence 
the correlation magnetic field of the fiuctuations in A/2 
contains information about spin relaxation, 

BcAh 0o/^' = <lyo/DT, , (6) 

with Is ~ \J being the spin-relaxation length. From 
Eqs. (^) and (^), we now calculate the spin-relaxation 
time 

~ (B,,G/S,,A/,)/r . (7) 

All parameters in this expression have been measured in 
our experiment. In particular, it does not depend on the 
diffusion constant. From the measured G(/?) and A/2(/?) 
data we find B^q ^ 0.01 T and B^a/^ ^ 0.03 T (in 
analogy to Ref. 0), and from the width of the conduc- 
tance peaks in Fig. ^(b) we deduce an energy resolution 
of fi,r « ea X 80 /iV ~ 40 /ieV. Using these values, we 
obtain a spin-relaxation time of Ts '--^ 5 ps. It is much 
longer than the elastic scattering time r = 0.1-0.2 ps, 
which we determined according to Ref. |p^ . 

Another interesting spin effect may result from the 
presence of magnetic impurities in the emitter. The ex- 
change field of fully spin-polarized impurities leads to 
different disorder potentials for spin-up and spin-down 



electrons and thus infiuences the LDOS (T9j. As there 
may only be residual magnetic impurities, their density 
will be by orders of magnitude lower than the Si con- 
centration. Therefore, the spin-dependent features of 
the disorder potential vary on longer length scales than 
the spin-independent features. Slow variations in space 
correspond to small features on short field scales in the 
LDOS. Thus a small amount of magnetic impurities could 
explain why small features in the LDOS exhibit an irreg- 
ular spin dependence, while large features simply show a 
regular Zeeman splitting. To clarify the influence of mag- 
netic impurities on the LDOS, we suggest to deliberately 
dope the GaAs contacts of a double-barrier heterostruc- 
ture with magnetic impurities such as Mn | p0| . 

In summary, we used a spin-split impurity level as spec- 
trometer to study the interplay of Zeeman splitting and 
disorder in Si doped GaAs via resonant tunneling. Both 
the spin-up and spin-down components of the LDOS 
show fan-type resonances as a consequence of Landau 
quantization in the bulk semiconductor. While the main 
resonances in the LDOS exhibit a Zeeman shift that can 
be best described using a g-factor of 5 = —0.35, the 
smaller resonances show a complicated irregular behav- 
ior. Spin relaxation due to spin-orbit coupling is one 
possible origin of this deviation from Zeeman splitting. 

We gratefully acknowledge A. Forster and H. Liith 
for growing the double-barrier heterostructure and V.I. 
Fal'ko for useful discussions. 
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FIG. 1. (a) Current-voltage curves I{V) measured at 
B = 8.0 T (solid) and 14.2 T (dashed) at T = 20 mK. (b) 
Differential conductance G — dl/dV numerically obtained 
from the I{V) curves. The sketch shows the spin-resolved lo- 
cal density of states (left) in the emitter of a double-barrier 
heterostructure (shown to the right). 

FIG. 2. Spin splitting AVs = Vs^ — Vs^ of the spectrometer 
vs magnetic field. 

FIG. 3. (Color), (a) Color map of the current vs bias volt- 
age (step 7 /iV) and magnetic field B \\ I (step 10 mT) for 
T = 20 mK (white, / < 0.0 nA; black, / > 0.5 nA). The raw 
data has been shifted in voltage such that the first current 
step defines origin of the voltage scale. The image is dom- 
inated by fan- type resonances Ve- (b) Current vs magnetic 
field at Vi (dashed, 7i) and V2 (solid, h). 

FIG. 4. Reconstruction of I2 from Ji: Comparison of the 
measured current I2 (solid) with the reconstructed current /I* 
(dashed). The solid line at the bottom shows the difference 
A/2 = /2 — /I* as a function of magnetic field. The sketch 
shows the local density of states (left) and the spectrometer 
position for Vi and V2 (right) to illustrate the reconstruction 
principle. Annotations are described in the text. 

FIG. 5. Root-mean-square deviation o-a/2 of 

A/2 — I2 ~ sls a, function of the (/-factor in the emitter. 
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Figure 2, Koenig et al. 




Figure 3, Koenig et al. 
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Figure 4, Koenig et al. 
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Figure 5, Koenig et al. 



